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Abstract 

We consider a generation and an amplification of THz radiation in semiconductor superlattices 
under the action of microwave pump field. Electrons belonging to a single miniband of the super- 
lattice interact quasistatically with the pump field and dynamically with a signal THz field. Within 
this semiquasistic approach we derive elegant difference formulas describing absorption (gain) of 
the weak THz signal. We present an instructive geometric interpretation of the absorption formu- 
las which allows a search of optimum conditions for the gain employing only a simple qualitative 
analysis. Our theoretical findings contribute to the development of sources and detectors of THz 
radiation that are using nonlinear electric properties of semiconductor superlattices. 
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I. INTRODUCTION 

Semiconductor superlattices (SSLs) 1] have attracted growing attention in view of their 
unique electronic properties, which can be used for generation, amplification and detection 
of a high-frequency electromagnetic radiation 

una 

. A dc-biased SSL, operating in the 
conditions of single miniband transport regime and demonstrating the static negative differ- 
ential conductance, can be potentially used as an active element of THz field generator . 
However, the static negative differential conductance makes SSL unstable against a forma- 
tion of hi g h- fi eld electric do matas QH. The electric do mainS are beiieved to b c destructive 



for the THz gain in SSLs 



. Currently the main focus is on the possibilities to overcome 



this drawback within the scheme of dc-biased SSL , However, schemes 

of THz superlattice devices with ac pump fields are also under discussion [3, [lj, 15, 
A strong microwave field E p (t) = E\ co${u)\i) pumps the SSL and a desirable signal field 
E s (t) = E 2 cos(u 2 t) has a higher frequency oj 2 > iO\- Because for typical SSLs the charac- 
teristic scattering time r at room temperature is of the order of 100 fs, an interaction of 
the microwave fields with the miniband electrons is quasistatic {oj\t <C 1). Importantly, 
we have showed recently that such quasistatic pump field can completely suppress domains 
m SSLs Two distinct possibilities exist for the signal field: It can be also qua- 

sistatic (uj 2 t <C 1) or it cannot be described within the quasistatic approach if uj 2 t > 1. 
The later situation, i.e. uj\t <C 1 but uj 2 t > 1, can be called semiquasistatic interaction. 
Thus, the semiquasistatic approach is introduced to describe an amplification of THz field 
in SSL under the action of microwave field. We should note that while in the experiments 



113, 



14j the interaction of electrons with both pump and signal fields is still quasistatic and, 



moreover, it is probably domains-mediated, the ultimate goal of this activity is to reach a 
parametric THz generation without electric domains |lj|. Therefore, we can say that the 
semiquasistatic approach describes the real experimental situations as well. 

In this paper, we derive several beautiful formulas describing a high-frequency absorption 
(gain) of small-signal field in the presence of ac pump field, when the interaction of the 
fields with the miniband electrons in the superlattice is semiquasistatic. Along with obvious 
potential applications in the calculations of gain in superlattice devices, our theoretical 
approach provides a powerful tool for finding the correspondence between quasistatic and 
dynamic regimes in ac-driven superlattices. 
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Our main findings are following. We consider the response of electrons, belonging to a 
single miniband of SSL, to an action of the total electric field 



E(t) = E + Ei cos(cjit) + E 2 cos(u 2 t), (1) 

where E is the dc bias and E s = E 2 cos(cu 2 t) is the weak signal (probe) field. We should note 
that in the real devices E s is a mode of the resonator tuned to a desirable THz frequency. 
Here three distinct cases should be considered: (i) u 2 and u)\ are incommensurable, (ii) 
uo 2 = muoi/2 (m is an odd number) and (iii) u 2 = muj\ (m is an integer). The later case 
of the parametric cascading is most interesting physically Q| and we mainly focus on it in 
this paper. We define the dimensionless absorption of the weak ac probe field in SSL as 

T 

A(lu 2 ) = 2(V(t) cos(uj 2 t)) t = ^ J V(t) cos(cu 2 t) dt, (2) 



where V(t) is the scaled miniband electron velocity defined in the units of the maximal 

nn 

miniband velocity Vq |lll2j|, averaging (. . .) t is performed over time and T = 2nm/uu 2 = 2ir/ui 
is the common period for the probe and the pump ac fields. 

Starting from the exact formal solution of the Boltzmann transport equation, we represent 
the absorption A as the sum of three terms 

j^harm _|_ j^coh _|_ j^incoh 

Here A harm , A coh and A mcoh describe the absorption (gain) seeded by generation of harmon- 
ics, the parametric amplification of the probe field due to a coherent interaction of the pump 
and the probe fields, and the nonparametric absorption, correspondingly. 

The term A harm is just the expression for mth in-phase harmonic of the time- dependent 
current through SSL 

A harm = 2 + ^ ^fafy cos ( muJl t)} t , (4) 

where Ud c = eLE is the dc voltage, L = Nd is the length of SSL (d is the period of SSL 
and N is the number spatial periods), U ac = eLEi is the amplitude of ac voltage created by 
the pump field across SSL, the current I(t) is normalized to the maximal current in SSL, I , 
corresponding to the maximal miniband velocity Vo- The explicit expressions for the mth 
harmonic of the current, that determine the dependence A harm (Ud c ,U ac ), are well-known 
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jrl 18 1 for the case of a small pump amplitude, U ac <C 1, and Q for the case 



of an arbitrary £/ ac ). Note that A harm does not depend on the probe; it gives the main 
contribution to the absorption of a weak probe [3, Q • 

The expression for A harm is not specific for the semiquasistic limit because of its indepen- 
dence on uj 2 . However, our main finding is that the absorption components A coh and A incoh 
can be represented within the semiquasistatic approach using the specific terms of quantum 
derivatives as 

, . / ' I ET (U dc + Uac cos(wit) + Nhw 2 /e) - I ET (U dc + U ac cos(c^) - Nhio 2 /e) 
A = eU 2 (^ cos(2mc 1 t) 

(5) 

I ET {U dc + U ac cos(o;i*) + Nhuj 2 /e) - I ET (U dc + U ac cos(c^) - NHuj 2 /e) ' 



A incoh = ^ 



2Nhu 2 



where U 2 is the amplitude of small-signal voltage and 

1 ] i + (u/u c y 1 ] 

is the Esaki-Tsu voltage-current (UI) characteristic (U c = HN/er is the critical voltage and 
I ET is normalized to the maximal current Iq = 2I ET (U = U c ) oc Vo). 

Importantly, following Eq. (0) in order to find the incoherent absorption in SSL at arbi- 
trary high frequency u 2 we need to know only 

I dc (U dc ) = (l ET {U dc + U ac cos(uj 1 t))) t , (8) 

that is, the time-averaged current induced by the quasistatic field (voltage). For a given 
amplitude of the ac voltage U ac , the dc current I dc is a function of only dc bias U dc . It is 
easy to calculate or to measure the modifications of UI characteristic caused by the action 
of microwave (quasistatic) field I21I 

On the other hand, finding of the coherent component of absorption at the high frequency, 
which corresponds to the mth harmonic of the pump frequency (u 2 = mui), is reduced 
following Eq. (JHJ) to the calculation of the 2mth harmonic of the current within the quasistatic 
approach. This result is especially important, because the term A coh can be responsible for 
the parametric gain in ac-driven SSLs [16]. Thus, one needs to know 

V^ rm (U dc ) = 2 (l ET (U dc + U ac cos^t)) cos(m^)) t . (9) 



(Note that in comparison with the expression (0J) we have in Eq. Q: I(U) — *> I ET (U)). For a 
given pump amplitude U ac the amplitudes of harmonics of the quasistatic current V^ rra are 
functions of only dc bias Ud c - We underline that quasistatic calculations of harmonics V^ rm 
are among the standard theoretical tools to explain experimental results on the frequency 
multiplications in the microwave frequency band [23, 24]. 

In this paper we will present simple but instructive geometric interpretations of the 
formulas (0) and © allowing to find the high-frequency absorption components directly 
using only the ruler and the knowledge of the quasistatic curves V£ arm (Udc) and Idc{U^}. 
These geometric interpretations provide a powerful tool in a search of the optimum conditions 
for the gain which employs only a simple qualitative analysis. 

If the frequency of the probe u>2 incommensurates with the frequency of the pump uii, 
the expression for the absorption at 002 consists of only one term - the incoherent absorption 
A = A mcoh with semiquasistatic limit ©• Additionally, we will also show that within the 
semiquasistatic approach the difference formulas similar to Eqs. ©-© can be written for 
the absorption at half-harmonics 0J2 = mui/2. 

The organization of this paper is as follows. We start with the calculations of absorption 
for arbitrary pump and probe frequencies solving the Boltzmann transport equation in 
sec. HU The derivation of the absorption formulas (j3J) and (jHJ) within the semiquasistatic 
limit is presented in sec. IIHI In the sec. IIV1 an additional limitation on the frequency of the 
probe is imposed and the quasistatic formulas for the absorption are obtained. Section Mis 
devoted to the presentation of useful geometric interpretations of the semiquasstatic results 
and their numerical verifications. The summary of the work, as well as a brief discussion 
of the place of our findings among other works devoted to the use of quantum derivatives 
in solid-state systems, are presented in the final section. Appendixes [X] and [B] are devoted 
respectively to the derivation of the semiquasistatic limit for the coherent component of 
absorption and to the derivation of the quasistatic limit for the harmonic component of 
absorption. 



II. ABSORPTION IN THE GENERAL CASE 

In this section we will represent the expressions for the components of absorption that 
follow from the exact formal solution of the transport Boltzmann equation with a constant 
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relaxation time j3, Using this exact solution, we found that for general case 

of commensurate frequencies uj\/uj2 = n/m (n,m are integers and n/m is an irreducible 
fraction) the total absorption A(u2), Eq. (J2J) has the form 



^0 t n^i -r i,2^2)i nm 



{i,ia=— do j"=— oo 

(Oq + ZiCUi + Z 2 Co> 2 )r 

1 + (fi + ^ + / 2 cu 2 ) 2 r 2 
where /3j = fli/uJi, = edEi/h (i = 1,2), and fio — edE /H is the Bloch frequency. To 
consider the absorption at mth harmonic of the pump frequency one should fix n = 1; for 
the case of absorption at half-harmonics n should be equal to 2. 

First, we consider the absorption at harmonics of the pump field. In the limit of the weak 
probe (/?2 C 1), we need to take only certain combinations of the Bessel functions indexes. 
As a result we represent the absorption as a sum of three terms (^tJ 

A = A harm + A coh + A incoh + 0(/3f), (11) 
' 1 + (iZ + <Co>i)^r z 

t=— oo 

ft y T/aUr / q \ j i a m (^ + /^i)r 

Z=— oo 



A- = f _£ M-a™(A) - WA)1 t/^iv (13) 

(fi + ^1 + ^ 2 )t (O + /to"i - uj 2 )t 



A incoh = Jf(j3 



2 ' ^ [l + (fi + Zu; a + cu 2 ) 2 r 2 1 + (fi + Z Wl - w 2 ) 2 r 2 

l=—oo 



(14) 



Second, for the absorption at half-harmonics of the pump frequency we get from (JTUJl the 
expression for A(u> 2 ) as the sum of two terms 

A h = A c h oh + A incoh + 0{(3l), (15) 

where 

x? = § ± m) k-<a) - j«-(a)1 ; («o 

!=— oo 

and A mcoh is given by the formula (|14|) . 

Finally, it can be shown that if u 2 /uji is some irrational number, then the absorption is 
completely determined by the incoherent component: A(u 2 ) = A mcoh . 

Now we turn to the consideration of semiquasistatic limits of different absorption com- 
ponents. 



III. DERIVATION OF THE SEMIQUASISTATIC FORMULAS 

The derivation of semiquasistatic formulas is based on the use of the following integral 
representation of the Esaki-Tsu characteristic yj 

oo 

/ET( "» = Tvk? s \ / exp (-7) sk ^ t)dL (17) 





28|. 



and the asymptotic saddle-point method 

We start with the consideration of semiquasistatic limit for the term A mcoh , Eq. ffT^j) . We 
rewrite A incoh in terms of I ET JT^) as 



30 



A incoh = yH J i^i) V ET i^o + iwi + u 2 ) - I ET (Q + lu, - co 2 )] 

l=—oo 



00 

00 



^ J exp f-^\ jr J?(Pi){ s in [(^o + lui + ^2)*] - sin [(Q + lw x - uj 2 )t] } dt. 

Q i=— OO 



i=— OO 

Using the addition formula 



sin [(fi + i^i + oj 2 )t\ = sin [($7 + L0 2 )t] cos(lu>it) + cos [(f2 + L0 2 )t\ sm(lujit) 

and taking into account that the sum with sin(/a; 1 t) is equal to zero, we get 

00 

Al ncoh = fh f (J_\ J2 Jf(fr){ sin [(Oo + iwi + wa)t] 

00 

- sin [(ffo + lu)\ - to 2 )t] I = ^ y exp ^0 f 2^ sin 



x I sin [(fi„ + w 2 )t] - sin [(O - u 2 )t] } dt. (18) 
In Eq. (JT5J) we also have used the equality j^J 

£ Jf (ft) cos(Za;it) = Jo (2/?! sin ^ . 
As a next step, we substitute the integral representation for Jq(z) 



Jq{z) = — I cos(z cos 6) d0 



7T 





in ([T%]) and get 



.4 



incoh 



A 

2vrr 



cxp 



7T 

J cos ^2/5i sin cos d6 j sin [(fi + co> 2 )t] 



o 

oo 



- sin [(Q - cj 2 )*] } dt = 7^- / exp 



sm 



^ sin(u; 1 t/2) ^ „ . 

fil C0S6 + ^ + ^2 t 

Uit/2 



sm 



„ sin(u;it/2) _ . 
fii — v , 7 cos6 + fi -o; 2 t 



uit/2 



dQ dt 



W1T<1 



A 

47TT 



2tt 



y exp ^— — ^ y | sin [(fix cos 6 + Q + u 2 ) t] — sin [(fix cos 6 + Q — uj 2 ) t) j dO dt. 



(19) 

In we first used simple formula 2 cos a sin /3 = sin(o; + /3) — sin(a; — 0) and then took 
into account that for uj-\t <C 1 the main contribution to the integrals with the exponential 
factor exp(— t/r) = exp(— c^x^/^i 7- ) comes from the terms formally satisfying u)\t — ► (the 
saddle-point method). Finally, using Eq. ()17j) we make the inverse transformation to the 
Esaki-Tsu current and obtain 



.4 



incoh 



fa 

4:71 



2- 



y [I st (fix cos e + n + u 2 ) - i ET (fix cos e + n - u 2 )] dQ. 



(20) 



Now it is easy to see that after making the transformation from frequency to voltage variables 
and using the Esaki-Tsu characteristic in the form (|7jl. we get for A m the formula (JBJ. 

Applying the same method, one can obtain (see Appendix EJ in the semiquasistatic limit 
the following expression for the coherent component of absorption 

2tt 

Xcoh fa 



A L 



2l 
4tt 



J [/ ST (fii cos e + n + uo 2 ) - i ET {n x cos e + n - u 2 )] cos(2me) t/e, (21) 



which is equivalent to the formula ©• 

We turn to the discussion of the semiquasistatic limit for absorption in the cases of 
incommensurable frequencies and half-harmonics. As we have found in the previous section, 
the absorption for incommensurable frequencies is just A mcoh defined by Eq. (|14|) . Therefore, 
its semiquasistatic limit is just the formula (J2Uj) . 

Next, in the general case of arbitrary values of uj\t and uj 2 t, the coherent component of 
absorption at half-harmonics A c h oh has the same form as the coherent component of absorp- 
tion at harmonics ([13)1 . if one makes a formal change 2m — ► m (cf. Eqs. ([IB}) and ([13])) . 
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Therefore, the quasistatic limit of A c £ h coincides with Eq. (jSJ) after the substitution 2m — > m. 
However, we should note that the consistent proof of this statement is not trivial because 
the derivations of semiquasistatic formulas for odd and for even values m have different ways 
(in the case of half-harmonics m is the odd number but in the case of harmonics 2m is the 
even number). 



IV. QUASISTATIC LIMIT 

To find the quasistatic limit of the components of absorption we have to impose the 
additional limitation on u 2 : lo 2 <C t _1 . For A harm we obtain from Eq. (|12j) in the quasistatic 
limit 

A harm = I j jET ^ + ^ cog Q ) cos ( m @) d Q ( 2 2) 


Of course, this expression is identical to Eq. ©. Formally, Eq. (}2~2~|) follows from (j3J) after 

the substitution I{U) —> I ET {U), but the correct and consistent derivation of this result still 

requires some algebra (see Appendix |B|) . 

To get the quasistatic limit for A tncoh and A coh we need to consider the limit uo 2 t <C 1 

in the Eqs. ((201) and (EI}- In this limit the finite difference /^(f^ cos 6 + Q + uj 2 ) - 

I ET (fl 1 cos 6 + Q — uj 2 ) goes into the derivative 

I ET (ni cos 9 + Sl Q + lo 2 ) - J £T (^i cos 9 + Q - lo 2 ) gJ £T (^ 1 cos9 + ^q) 

2uj 2 dVlo 

and hence A coh and A mcoh become 

, n 2 2 fdi ET (n lC ose + Q )_,^^ Aincoh n 2 J di ET (n lC osQ + n ) 



A con = — / ^— : — ^ cos(2m6) dQ, A incoh = — / : — ^ dQ. 

2tt J dQ 2tt J dfl 



These formulas in terms of voltage and current, Eq. (jJJ), have the following form 

A coh = U 2 ^j- (I ET (U dc + U ac cos(^t)) cos(2mo; 1 *)> t , (23) 

A mcoh = U 2 ^j- (l ET (U dc + 11^ cos(u 1 t))) t . (24) 

Thus, in the quasistatic limit the incoherent and the coherent components of absorption at 
mth harmonic are proportional respectively to the derivative of dc component (JBJ) and to 
the derivative of 2mth harmonic (151 . 
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On the other hand, the formulas for the absorption components can be derived completely 
within the quasistatic approach. We define the absorption of quasistatic field (voltage) as 
A qst {ui 2 ) = {Iet{U) cos(uj 2 t)) t , where the Esaki-Tsu UI characteristic Iet{U) is given by 
Eq. (J7J). For a weak signal we expand the current in the Taylor series I E t{U) ~ Iet{U p ) + 
I' ET {U P ) x U s (here prime means the derivation with respect to U, U p = U dc + U ac cos(c<j 1 t), 
U s = U 2 cos(u 2 t)). Substituting this expansion to the definition of A gst we have 

A« st (u 2 ) = 2(I ET (U p ) cos(u 2 t)) t + (I' ET (U P ) cos(2u 2 t)) t U 2 + (I' ET (U p )) t U 2 (25) 

Making comparison of Eq. (|23j) with other formulas from this section, we see that in the limit 
uj 2 t — > the semiquasistatic formulas for the components of absorption correctly reproduce 
the quasistatic results. 



V. GEOMETRICAL INTERPRETATION AND NUMERICAL VERIFICATION 

To discuss the geometric meaning of incoherent and coherent components of absorption 
it is convenient to present A mcoh and A coh explicitly in terms of zero and even harmonics of 
the quasistatic current I ET . Combining (j5J) and (jUJ) with (jHJ) and (jSJ) we have 

A mcoh = ^ [I dc (U dc + Nhw 2 /e) - I dc (U dc - Nhuj 2 /e)} (26) 

ZCo>2 

A coh = ~± [V 2 t rm (U dc + Nhw 2 /e) - V 2 t rm (U dc - NfHu 2 /e)] (27) 

Alternatively, the expressions for the absorption components in terms of frequency (|2*U|) . (|2*Tj) 
can be used. 

As is evident from Fig. Q the geometric representation of the quantum derivative (finite 
difference) of I dc is the segment with the length defined by the frequency of probe field uj 2 . 
The ends of the segment belong to the curve I d c{U dc ) and their locations are determined by 
the location of working point. Slope of the segment determines A mcoh within the semiqua- 
sistatic approach. In the quasistatic limit, a finite difference becomes an usual derivative 
and therefore A %ncoh is determined just by the slope of the tangent line to I dc {U dc ) curve at 
the working point. In a similar manner, the coherent component of absorption A coh at the 
frequency uj 2 = mui is determined by the slope of the segment with the ends belonging to 
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FIG. 1: [Color online] Geometrical meaning of the incoherent component of absorption within semi- 
quasistatic and quasistatic approaches. Time-averaged current Idc under the action of quasistatic 
pump {<jJ\t = 0.1, U ac /U c = 0.2) vs dc voltage Udc- If we choose the working point at Udc = 0, then 
the dotted segment [red online] corresponds to the finite difference (quantum derivative) of a weak 
probe at L02 = lltoi and the dashed straight line [blue online] corresponds to the derivative. The 
slopes of these strait lines determines A mcoh respectively in semiquasistatic and quasistatic cases. 

the curve V^ rm (U dc ) (Fig.EJ). Respectively, the slope of the tangent line to V^ rm (U dc ) at 
the working point determines A coh in the quasistatic limit. 

Using these two geometrical pictures it is possible, in principle, to optimize a location 
of the working point (i.e. the value of applied dc bias) in order to obtain the maximal 
small-signal gain in SSL in the conditions of suppressed domains, when both m = lOijuJx 
and the pump amplitude U ac are given. However, such an analysis goes beyond the scope of 
the present paper. In the rest of this section we will focus on the case of dc-unbaised SSL 
{Udc = 0), which corresponds to the choice of working point at the origin of coordinates in 
Figs. [U and El Such a choice guarantees the suppression of electric domains |l^ : a slope 
of the dependence Id c on Ud c is positive practically for all U ac in the limit Udc (for 
examples, see Figs. [H HI EJ). 

To illustrate the usefulness of the geometric interpretation in finding of schemes allowing 
gain in the conditions of suppressed domains, we will consider two characteristic examples 
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FIG. 2: [color online] Geometrical meaning of the coherent component of absorption A coh within 
semiquasistatic and quasistatic approaches. The 2mth harmonic of quasistatic current V 2 h 1 ^ rm ' vs 
dc voltage Udc for m = 7. The current is induced by the quasistatic pump (uj\t = 0.14) with 
the amplitude U ac /U c = 8. Slope of dotted segment [red online] determines the (negative) finite 
difference for a weak probe with 002 = 7oo\. Slope of dashed straight line [blue online] determines 
the derivative. Working point is chosen at Udc = 0. 

(earlier the solutions of these problems have been announced without proofs in [16]). 

Problem I. Suppose that 002/ is an irrational number and therefore only the incoherent 
absorption A mcoh contributes to the gain in SSL. Employing the geometrical interpretation 
of A m , we immediately see that if the slope of the tangent line to Id c curve is positive, then 
the quantum derivative is also positive for an arbitrary u>2 (see, e. g., Figs. Il|4|) . In physical 
terms, this means that a small-signal gain in the unbiased SSL is impossible in the conditions 
of suppressed domains, if the frequencies of pump and probe are incommensurable. 

Problem II. Consider the absorption at a probe frequency that is a half-harmonic of 
the pump, U2/0J1 = m/2. The total absorption is the sum of the incoherent component 
A mcoh and the coherent absorption at half-harmonics A c £ h (see Eq. [Tty . The geometric 
interpretation of A c £ h is similar to the interpretation of A coh : One needs only to consider 
yharm ag a f unc ^i on Q f jj ^ instead of a dependence V 2 h ^" m {Udc) ■ Obviously, the dependence 
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FIG. 3: [color online] Geometric interpretation for the coherent component of absorption at half- 
harmonics A c £ h within semiquasistatic (dotted [red online] segment) and quasistatic (dashed [blue 
online] line). The curve represents rath harmonic of quasistatic current Vj^ arm vs dc voltage Udc 
for u>2 = 7u>x/2, uj\t = 0.14 and Oir = 8. The zero slopes of the straight lines show that Af^ h 
calculated for Udc = is zero. 

of any harmonic of the quasistatic current on the applied dc bias, i. e. V^ arm (f/ (ic ), is always 
symmetric about the vertical line Udc — Therefore, the choice of the working point at 
Udc = results in zero value of both the derivative and the quantum derivative in this point 
(Fig. El). Thus, in the case of dc-unbiased SSL A c ^ h is always zero. On the other hand, 
as we have saw earlier, A mcoh is always positive in the conditions of suppressed domains. 
Summing up we can conclude that a small-signal gain in unbiased SSL is impossible in the 
conditions of suppressed domains, if the frequency of probe is a half-harmonic of the pump 
frequency. 

These two conclusive findings allow to narrow the class of possible schemes for the su- 
perlattice THz emitters with ac-pump, which can operate in the domainless regimes. The 
proper scheme should exploit a negative absorption at some particular harmonics of the 
probe, when negative values of harmonic (A harm ) and coherent contributions to the net 

n 

gain can overcome always positive incoherent absorption component |16j . Therefore, in our 
analysis we return to the case co^/^i =integer and compare A incoh and A coh in conditions of 
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FIG. 4: [Color online] Incoherent component of absorption for different strengths of ac pump. 
Both difference and usual derivatives are simultaneously small and approach each other for the 
strong pump (U ac /U c = 2) [blue online], in contrast to the case of weak pump (U ac /U c = 0.2) [red 
online] where they are quite different. Here uj\t = 0.14 and m = 11. 

quasistatic, semiquasistatic and dynamic approaches. 

Semiquasistatic and quasistatic results have a visible difference for higher harmonics 
as is evident from Figs, d and El (here the quasistatic condition uo 2 t <d 1 is not satisfied 
for the probe field). However, in the case of low-order harmonics the difference between 
semiquasistatic and quasistatic results is very small. This fact is illustrated in Fig. for 
A coh , where the difference between the derivative (tangent line) and the quantum derivative 
is very small for m — 3, but it visible increases for m = 7. 

The dependencies of A incoh and A coh on the amplitude of pump field, calculated using 
semiquasistatic (solid) and quasistatic (dashed) formulas for large m, are shown in Figs. El 
and [3 Here A mcoh is always positive and its value decreases with an increase of the pump 
amplitude. Fig. El demonstrates that the quasistatic calculations significantly overestimate 
A mcoh for a small pump U ac /U c < 1 (see also Fig. [TJ. However, the agreement between 
semiquasistatic and quasistatic approaches in finding A mcoh becomes good for a strong pump 
(Fig. IHJ). Such paradoxical behavior is explained by the effect of shift of the maximum of UI 
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FIG. 5: [Color online] Comparison of the coherent absorptions for different values of m = o^/t^i- 
Quantum derivative [blue online] approaches normal derivative (the tangent line, not shown) for 
m = 3; they are quite different for m = 7 [red online]. Quasistatic pump with lot = 0.14 has the 
amplitude U ac /U c = 8. 



characteristic under the action of quasistatic ac field [30J, |3l| (Note: in our case the pump 
is still quasistatic). As we see in Fig. HJ for a large pump UI characteristic is rather closed 
to the linear in the interval of dc voltages that is utilized for the calculation of quantum 
derivative, what is in a sharp contrast to the case of large ac pump where UI characteristic 
demonstrates a strongly nonlinear behavior in the same interval. Thus, for a strong pump the 
finite difference is close to the tangent line, providing a small difference between the results 
of semiquastatic and quasistatic approaches for A mcoh . As for A coh , the difference between 
semiquastatic and quasistatic calculations is quite visible for large m (Fig. EJ). Moreover, 
the quasistatic approach not only overestimate gain in this case, it can give even a wrong 
sign of A coh in some interval of ac pump amplitudes (Fig. |7J). 

We also calculated numerically A mcoh and A mcoh using formulas ()13p. (fT4"|) valid for arbi- 
trary values of uj\t and 0J2T. Results of these calculations in comparison with semiquasistatic 
calculations are shown in Figs. |HJ El for the case of moderate pump strength and in Figs. lTUllTTI 
for a strong ac pump. For all calculations, including those which are not shown in these 
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FIG. 6: Dependence of A m on the amplitude of ac pump calculated within semiquasistatic 
(solid) and quasistatic (dashed) approaches for m = 9 and oj\t = 0.14. 
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FIG. 8: Dependence of A m on oj\t calculated within semiquasistatic (solid), dynamic (dotted) 
and quasistatic (dashed) approaches for m = 5 and pump amplitude £l\T = 2.0. 

figures, we always found an excellent agreement between results of dynamic (dots) and 
semiquasistatic (solid) approaches for small uj\t. Surprisingly, the semiquasistatic formulas 
work quite well even beyond their formal range of validity, i.e. not only for small values 
of U\T. For the moderate pump, the semiquasistatic approach works well up to uj\t ~ 2 
(Figures |Hl El)- Even for a strong enough pump, the semiquasistatic approach still gives a 
correct smoothed-out curve, while it does not exactly reproduce small oscillations in the 
dependence of absorption on the pump frequency (Figures llOll 1(1 . 

VI. CONCLUSION 

In summary, we calculated a small-signal absorption (gain) of a high-frequency field in 
semiconductor superlattice driven by a quasistatic ac field. In particular, our theoreti- 
cal approach well describes practically interesting cases of generation and amplification of 
THz signal in superlattice devices pumped by a microwave field. We presented transpar- 
ent geometrical interpretations for the coherent and the incoherent contributions to the net 
absorption at the different ratio of the pump (ux) and the signal (0*2) frequencies. We 
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FIG. 9: Dependence of A coh on uj\t calculated within semiquasistatic (solid), dynamic (dotted) 
and quasistatic (dashed) approaches for m = 5 and Q\t = 2.0. 
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FIG. 10: Same as in Fig. H but for = 5.1. 
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FIG. 11: Same as in Fig. El but for = 5.1. 



demonstrated that the use of these geometric pictures sufficiently simplifies analysis: It is 
easy to find gain for a signal with an arbitrary u 2 using the electric characteristics of su- 
per lattice formed under microwave fields. In particular, we proved that for iOiju)\ being 
some irrational number a small-signal gain at io<i is possible only in the condition of static 
negative differential conductance. Next, we also found that an absorption at an arbitrary 
co>2, that equals to some half-integer of the pump frequency, is always zero in an unbiased 
superlattice. At least but not last, we presented the universal analytic procedure for find- 
ing the behavior of different physically important variables in the quasistatic limit starting 
from the exact solution of the Boltzmann transport equations for superlattices. Numerical 
calculations confirmed our analytic findings and also demonstrated that the semiquasistatic 
formulas for absorption work reasonable well even beyond their formal range of validity. Our 
results can be generalized to other physical systems where the miniband transport regime 
exist, such as carbon nanotubes or dissipative optical lattices 

Although in this work we used only a single scattering time, our results obtained within 
semiquasistatic approximation can be also applied to a more realistic case of two different 
relaxation times for electron energy and electron momentum. Really, within the approxi- 
mation of two relaxation times the static voltage-current characteristic of superlattice still 
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has the Esaki-Tsu form |34(. Because in the semiquasistatic approach only the (quasi)static 
volt age- current characteristic plays a role, we speculate that a variation of the ratio of re- 
laxation times can lead only to a some change in absolute value of the absorption but not 
in its dependence on the amplitude of pump field. Preliminary numerical simulations of the 
superlattice balance equations 3, 3] confirm these speculations 3^|. The detailed consid- 
eration of THz absorption with the use of two different relaxation times will be published 
elsewhere. 

It is instructive to consider our activity from the viewpoint of its place among other 
works devoted to the use of the quantum derivatives or, other words, to the use of the finite 
differences in expressions for currents. In the theory of semiconductor superlattices, the 
quantum derivative has been first introduced in [36] for calculations of the absorption of 
weak monochromatic ac field in dc-biased superlattice. Later the formula has been used in 
the calculations of the response function of THz superlattice detector |37( and in the fin ding 



of THz gain in the models of generalized Bloch oscillator with suppressed domains 
Interestingly, these formulas for superlattices appears to be same as the formulas describing 
the response of two-terminal structures to the action of dc and ac voltages 40]. The latter 
are widely used for calculations of signal mixing in the Josephson junctions operating in the 
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qusiparticle transport regime and in other single-barrier tunnel devices 

To the best of our knowledge this paper is the first work, where the quantum deriva- 
tives naturally appeared in a description of microstructure's response to a bichromatic field. 
It is an exciting problem to understand whether our semiquasistatic approach, developed 
for superlattices operating in the miniband transport regime, can be generalized to other 
tunnelling structures. 
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APPENDIX A: SEMIQUASISTATIC LIMIT FOR A coh 



In this Appendix we consider the semiquasistatic limit for A coh . We rewrite A coh ()13|) in 
terms of I ET (I17|) as 

A coh = t^Y, Jl(/3l)[Jl^ m Wl)-Jl + 2 m Wl)]I ET ^0 + ^l) 
l=—oo 
00 

= ^ J e*p(~) Ji(Pi)[Ji-2m(Pi)-Ji + 2UPi)]s™[(n + lui)t]dt. (Al) 

q l = — 00 

After changing indexes, Z — > — Z, in the first term of Eq. (jAljl we obtain 

00 

A00 fc = ^ /" exp^-^J ]T J/(A)^+2m(/3i){sin[(n -/a;i)t] - sin [(n + L^t] } dt. 

Q i = — OO 

(A2) 

Using the addition formula for sines we get 

00 

A co?i = /" eotpf~) JK/5i)^+2 m (/3i)sin(Z^t)cos(fiot)^. 

i=— CO 



(A3) 



To make the summation over Z we use the identity 

00 / 
2^ sm(la)Ji(x)Ji +m (y) = sin | marcsin 



xsma 



l=—oc 



\/ x 2 + y 2 — 2xy cos a 
x J m ^a/x 2 + y 2 — 2xy cos aj (A4) 

and get 

00 

^cofc, _ _$2l J eX p ^ \ s i n ^2marcsincos -^-^j ^2(3i sin cos(fioO 



(A5) 

Now we can take into account that for uj\T <C 1 the main contribution to the integrals with 
an exponential factor comes from the terms formally satisfying io\t — > 0. In this case, we 
can use the estimations 

h m \ 2fh sin - Jim (fii*) , 
21 



because 



44] 



sin ( 2m arcsin cos ) ~ (— l) m si^mu^t) = (— l) m sin^t), 



7T . 

arcsin(l — z) = — — v2z, for \z\ <C 1. 



Substituting these estimations to (|A5|) we derive 

oo 

A coh = _izlLA y exp f-^j sin(a; 2 t)cos(^)^2m(^i^)^ 



o 

oo 



^ 2 y exp ^— "J J 2m(^i*) [sm(u 2 t - Q t) + sin(cj 2 t + Q t)] dt. (A6) 
o 



2r 

Using the integral representation for J m (z) 

IT 

1 ir 

Jm{z) = — / cos(2; sin — mQ) dO, (A7) 

y 



7T OO 

■l) m „ / / ft 



we have 

A coh = ^—p 2 J J exp f~J cos(fiitsine-2m6) 
o o 

x I sin [(O - w 2 )t] - sin [(O + w 2 )t] } dt d0. ( A8) 

In this expression we use the formula for a product of sines and cosines and as a result we 
get the sum of eight integrals. These integrals are 

IT IT 

y cos [Oitsine±(O -cj 2 )t] sin(2m9) dQ = J cos [Sl l tsinQ±(Slr ) +u 2 )t\ sin(2m6) dB = 0. 



o 



and 



7T 

J | sin [O x t sin 6 + (fi =F w 2 )t] - sin [f^t sin 6 - (Sl =F u 2 )t] j cos(2m6) d6 
o 

= 2(-l) m y sin [Sltf cos 6 + (fi =F u 2 )t] cos(2m6) d6. (A9) 
o 

Now ()A8|) takes the form 

7T OO 

A co/l = -J—P2 J J exp f — "J { sin cos © + (^0 - wa)*] - sin [fijt cos 9 + (fi„ + Wa)*] } 





x cos(2mQ)dtdQ, (A10) 
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which finally allow us to make the inverse transformation by means of (|17|) to the Esaki-Tsu 
current as 

A coh = — fo [ [J BT (fii cos 6 + tt + u 2 ) ~I ET {^i cos6 + tt -u 2 )] 
4n J 
o 

x cos(2m6)de = — (3 2 I \l ET (^li cos(cjii) + Q + uj 2 ) - I ET (Qi coa(uit) + Q -uj 2 
An J 
o 

x cos(2muJit)dt. (All 
That is the desirable semiquasistatic representation for A coh . 

APPENDIX B: QUASISTATIC LIMIT FOR A harm 

Here we consider the semiquasistatic limit for A harm . It is convenient for our purposes to 
change indexes I — > — I in the first term of (|T2"j) and rewrite A harm in the following form 



oo 



A harm = ^ J^fc) J^fc 



l=—oo 



(Qq — luj\)T (Oq + Ioj\)t 



1 + (fi„ - L;i) 2 r 2 l + (fi + ^i) 2 r 2 . 



(Bl) 

Now we need to rewrite A harm in terms of I st ((T7j) as 

oo 

A harm = JiWJi+nW [(-i) m i ET (n - luji) + i ET (n + lull 



l=—oo 

oo 

If {I 

T 



J exp f-^j jr JK/5i)^ +m (/3i){(-l) m sin [(Q - lujt] + sin [(n + lujt] } dt 

Z=— oo 

oo 

= i /" expf-~J ]T J,^) J J+m (A){[(-l) m + 1] sin(fiot) cos(/ Wl t) 

Q ^ i = — OO 

+ [(-l) m -l]cos(O i)sin(io;it)}^ (B2) 
We can take the sum over I in (jB2|) using Formula (jA4|) and the identity 2^] 

oo / 

2 cos(la)Ji(x)Ji +m (y) = cos m arcsin = 
i= _ oo ' y a/x 2 + y 2 - 2xycosa^ 

x J m [\/x 2 + y 2 - 2xy cos aj . (B3) 
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As a result we have 

oo 



A harm = i J exp ^-^J { [(-l) m + 1] sin(Oot) cos (m arcsin cos J m ^2/3i sin 
o 

+ [(-l) m - 1] cos(fio^) sin arcsin cos J m f 2/3 x sin } eft. (B4) 

For cJiT <C 1 the main contribution to the integrals with an exponential factor comes from 

the terms formally satisfying u\t — > 0. Therefore, the following limits 

. / . Uit\ . rmx 
sm m arcsin cos ~ sin , 

/ Wit\ mm 

cos I m arcsin cos J ~ cos (B5) 

J m ^2/3i sin ~ J m (ftji) 
are important. Substituting (jB5|) in ()B4|) . we obtain 

oo 

A harm = l -J exp(~) {[(-l) m + l]sin(fi t)cos^J m (fi 1 t) 



+ [(-l) m - 1] cos((V) sin J m (£lib)j dt. (B6) 

Now odd and even values of m must be considered separately. 

If m — 2k + 1 (k = 0, 1, 2, . . .) is an odd number, then Eq. (|B6|) becomes 



oo 

A harm = ^—^ — / expl--Jooa(Q t)J m (n 1 t)dt 



o 

7T OO 



2( — l) k+1 [[ft 

exp I — ) cos(fio^) cos(f2]t sin O — mQ) dt dQ, (B7) 




itt j j V r 





where we also used the integral representation (|A7|) for J m . 

To simplify the expression (jB7|) we will need the following formula 



j cos(fV) cos(fiit sin 6 - mO) dB = (-l) fe+1 J sin[(fi + ^i cos Q)t] cos(m9) d6. 



(B8) 

To prove it we first convert the integral as 

7T 7T 

J cosing sin Q -mQ)de = J sin(fM sin 6) sin(m6) d0, (B9) 
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and then use the formula for a product of sines and cosines 

2 cos (Slot) sin sin 6) = sin[(fi - ^i sin Q)t] + sin[(fi + fii cos 0)t] . (BIO) 
Taking into account that 

IT IT 

J sin[(fi - fii sin 0)t] sin(m6) rf0 - j sin[(fi + fii sin 6)t] sin(mG) d6 
o o 

7T 

= 2(-l) k+1 J sin[(fi o + fi 1 cos0)t]cos(me)rf0, (Bll) 



we get (IB8D . 

Now substituting (jB8|) in ()B7|) and making the inverse transformation to the Esaki-Tsu 
current, we obtain 

7T OO 

A harm = A J j exp (--M sin [(O + fii cos6)t] cos(m6) dtdQ 



2vr 2tt/o;i 

= — I I ET (Q + fii cos 0) cos(m9) c/9 = — / I ET (Q + fii cos(wit)) coa(muit) dt. 

TIT J 71 J 



(B12) 



We turn now to the consideration of even harmonics. Combining (jB6|) for m = 2k 
(k = 1, 2, . . .) and the integral representation (jA7|) . we have 



k °° 

A harm = 2 ^~^ y exp f--^ sin(fi t) J ro (fiit) 



o 

7T OO 




IT J J \ T 





2(-lp f f f t 

exp I — I sin(f2 t) cos(f2it sin — mO) (it dO. (B13) 



Next stage of simplification includes several steps. First, we take into account that 

7T 7T 

J cos(fiitsin6 - mQ)dQ = j cos(^it sin 9) cos(rn9) d9. (B14) 



Second, we apply the formula for a product of sines and cosines. Finally we also use 

7T 

J | sin [(fi + fii sin 6)t] + sin [(fi - n x sin 6)i] | cos(m0) d6 
o 

= 2(-l) fc y sin [(O + fii cos 6)t] cos(m6) d9. (B15) 
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As a result, we can represent (jB13|) in the form suitable for the inverse transformation to 
the Esaki-Tsu current and obtain 

7T OO 

A harm = A J J exp sin [(O + fli cos6)t] cos(m6) dtdQ 



= \ J lET ( Q o + fti cos 0) cos(m6)rfe 



o 



— J I ET [VLq + fli cos(co>it)) cos(mciJit) (it 



The last equation coincides with the corresponding equation for odd harmonics. Thus, 
we conclude that in the quasistatic limit for both even and odd values of m 

2tt/uji 

A harm = ^ J J £T (fi + ^1 COS^t)) COs(mtU 1 t) dt. 


(B16) 
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